Search results

Two common fixed grid enthalpy methods used in the numerical modelling of phase change problems are the apparent heat capacity and the source based methods. In this paper, a general enthalpy method that includes as subsets both apparent heat capacity and source based methods, is derived. Following this, an optimal enthalpy scheme is identified. The superiority of the optimal scheme over the apparent heat capacity and the source based schemes is illustrated by solving sample phase change problems.

Cyclic phase change involves the successive freezing and melting of a region driven by a boundary temperature that cycles above and below the solid/liquid phase change temperature. In this paper, a recently proposed fixed grid phase change enthalpy method is modified and applied to cyclic solid/liquid phase change problems. The basic approach is demonstrated on application to a one‐dimensional, heat conduction controlled phase change. Then the method is used to investigate a cyclic phase change problem that involves fluid flow. The interaction of the melting and freezing with the phase change leads to some interesting predictions for the location and shape of the solid/liquid interface. The results also indicate that melting cycles are more effective than freezing cycles.

This paper aims to propose a semi-analytical benchmarking framework for enthalpy-based methods used in problems involving phase change with latent heat. The benchmark is based on a class of semi-analytical solutions of spatially symmetric Stefan problems in an arbitrary spatial dimension. Via a public repository this study provides a finite element numerical code based on the FEniCS computational platform, which can be used to test and compare any method of choice with the (semi-)analytical solutions. As a particular demonstration, this paper uses the benchmark to test several standard temperature-based implementations of the enthalpy method and assesses their accuracy and stability with respect to the discretization parameters.

The class of spatially symmetric semi-analytical self-similar solutions to the Stefan problem is found for an arbitrary spatial dimension, connecting some of the known results in a unified manner, while providing the solutions’ existence and uniqueness. For two chosen standard semi-implicit temperature-based enthalpy methods, the numerical error assessment of the implementations is carried out in the finite element formulation of the problem. This paper compares the numerical approximations to the semi-analytical solutions and analyzes the influence of discretization parameters, as well as their interdependence. This study also compares accuracy of these methods with other traditional approach based on time-explicit treatment of the effective heat capacity with and without iterative correction.

This study shows that the quantitative comparison between the semi-analytical and numerical solutions of the symmetric Stefan problems can serve as a robust tool for identifying the optimal values of discretization parameters, both in terms of accuracy and stability. Moreover, this study concludes that, from the performance point of view, both of the semi-implicit implementations studied are equivalent, for optimal choice of discretization parameters, they outperform the effective heat capacity method with iterative correction in terms of accuracy, but, by contrast, they lose stability for subcritical thickness of the mushy region.

The proposed benchmark provides a versatile, accessible test bed for computational methods approximating multidimensional phase change problems. The supplemented numerical code can be directly used to test any method of choice against the semi-analytical solutions.

While the solutions of the symmetric Stefan problems for individual spatial dimensions can be found scattered across the literature, the unifying perspective on their derivation presented here has, to the best of the authors’ knowledge, been missing. The unified formulation in a general dimension can be used for the systematic construction of well-posed, reliable and genuinely multidimensional benchmark experiments.

The axisymmetric steady‐state convective‐diffusive thermal field problem associated with direct‐chill, semi‐continuously cast billets has been solved using the dual reciprocity boundary element method. The solution is based on a formulation which incorporates the one‐phase physical model, Laplace equation fundamental solution weighting, and scaled augmented thin plate splines for transforming the domain integrals into a finite series of boundary integrals. Realistic non‐linear boundary conditions and temperature variation of all material properties are included. The solution is verified by comparison with the results of the classical finite volume method. Results for a 0.500[m] diameter Al 4.5 per cent Cu alloy billet at typical casting conditions are given.

Effective internal control over financial reporting (ICFR) should either prevent or enable correction of any material misstatement in a firm’s financial statements. Independent auditors, guided by professional standards, prepare an ICFR audit report, which provides a gauge by which the public can evaluate the reliability of a firm’s financial information. A firm manager may be tempted to misstate financial statements if he/she perceives a substantial reward for doing so. This study examines whether managers and firms are rewarded for misstating their financial statements in situations where there are incentives to do so, specifically, when an industry-leading peer is fraudulently inflating its reported earnings. The authors test to see if managers experience an increase in compensation as a result of misstatement. The authors also test to see if their firms benefit from misstating via changes to their cost of capital. Results suggest that neither managers nor their firms benefit from managing earnings by misstating financial statements. These findings are important, because a manager who ex ante understands that misstating will not lead to benefits personally or his/her firm is less likely to misstate in the first place.

The homographic approximation, in which the Heaviside step function is replaced by a continuous smooth curve, is applied to the enthalpy method for heat transfer problems with isothermal phase change. Both the finite difference and finite element implementations, based on the basic enthalpy, the apparent heat capacity and the source term formulations, are considered. A 1‐D Stefan problem of melting a solid is used as a test problem. The accuracy of the numerical solutions is measured globally using L2 error norms and comparison is made between the solutions using homographic approximation and those using linear approximation. The advantages of using homographic approximation are examined.

In the original effective heat capacity method the latent heat effect is approximated by a large effective heat capacity over a small temperature range. This method is simple in concept and easy to implement. However, it is so sensitive to the choice of the phase change temperature interval and the chosen time integration scheme that non‐convergence always occurs when it is implemented with implicit time integration schemes. Presents the cause and cure of non‐convergence in effective heat capacity methods and discusses sample results.

A control‐volume based method for the numerical calculation of axisymmetric incompressible fluid flow and heat transfer is presented. The proposed method extends the staggered grid approach to unstructured triangular meshes. The velocities are stored at the vertices and the edges of a triangle, pressure and temperature are stored at the vertices. Accordingly, velocities are interpolated in a quadratic way, pressure and temperature linearly. The accuracy of the proposed method is examined for a number of different testproblems. Compared to a linear interpolation scheme implemented in the same code, more accurate solutions and smaller computation times are obtained for the proposed quadratic scheme. The method was designed for and is about to be applied to the numerical simulation of crystal growth.

Proper selection of the stability parameter determines the accuracy of dendrite tip kinetics at a single crystal scale. Recently developed sophisticated phase field modelling of a single grain evolution provides evidence that this parameter is not constant during the process. Nevertheless, in the commonly used micro-macroscopic simulations of alloy solidification, it is a common practice to use a constant value of the stability parameter, resulting from the marginal stability theory. This paper aims to address the issue of how this inaccuracy in modelling crystal growth kinetics can influence numerically predicted zones of columnar and equiaxed dendrites and the macro-segregation formation.

Using the original authors’ micro-macroscopic computer simulation model of binary alloy solidification, the calculations have been performed for the Kurz-Giovanola-Trivedi (KGT) crystal growth kinetics with two different values of the stability parameter, and for two different compositions of Al-Cu alloys. The computational model is based on single domain-based formulation of transport equations, which are discretized on control-volume mesh. To identify zones of different grain structures, developing within the two-phase liquid-solid region, an envelope of columnar dendrite tips is tracked on a fixed non-orthogonal, triangular control volume grid. The models of porous and slurry media are used, along with the concept of the switching function, to account for diverse flow resistances in the columnar and equiaxed crystal zones. The numerical predictions are carefully studied to address the question of how the chosen stability parameter influences macroscopic structures of a cast, the most important issue from the engineering point of view.

The carried-out comprehensive numerical analysis shows that the value of the stability parameter of the KGT-constrained dendrite growth model does not have a direct significant impact on the macrosegregation formation. It, however, visibly influences the undercooling along the front, separating different dendritic structures and the size of the undercooled melt region where the equiaxed grains can develop. It also affects the amount of eutectic phase created.

To the best of the authors’ knowledge, this is the first attempt at estimating the influence of some inaccuracies, caused by possible ambiguities in choosing the stability constant of the KGT law, on numerically predicted macroscopic fields of solute concentration, the developing zones of columnar and equiaxed crystals and the macrosegregation patterns.

The purpose of this paper is to study thermal performance of metal foam/phase change materials composite under the influence of the enclosure aspect ratios (ratio of enclosure height: length). In this study, a compound metal foam/phase change material (PCM), which has been proved to be one of the most promising approaches for thermal conductivity promotion on PCMs, was used.

The PCM is considered initially at its melting temperature. The enclosure for all the cases has a constant volume with various aspect ratios. The left side of the enclosure is suddenly exposed to a thermal source having a constant heat flux, while the other three surfaces are kept thermally insulated. A two-dimensional numerical model considering the non-equilibrium thermal factor, non-Darcy effect and local natural convection was proposed. The coupling between velocity and pressure is solved using the SIMPLEC, and the Rhie and Chow interpolation is used to avoid the checker-board solutions for the pressure.

The effects of foam porosity and aspect ratio of the enclosure on the PCM’s melting time were investigated. The results indicated that enclosure aspect ratio plays a fundamental role in phase change of copper foam/PCM composites. For higher porosities, enclosures with bigger aspect ratios proved to led to optimal melting time. Besides, the best enclosure aspect ratio and foam porosity for a fixed-volume enclosure to have the shortest melting time are 2.1 and 91.66 per cent, respectively. However, for a specific amount of PCM inside a variable volume enclosure, the optimal melting time was for foam with ε = 95 per cent. The achieved results prove the great importance of selection of aspect ratio to benefit both conduction and convection heat transfer simultaneously.

The area of energy storage systems is original.